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Abstract An approximation for the inverse-Compton radiation power of electrons in the isotropic black-body 
photon field is presented. The approximation allows one to calculate inverse-Compton emissivity as integral over the 
energies of incident electrons rather than over the field photon energies. Such an approach allows for accurate modeling 
of IC emission of electrons with energy spectra being different from power-law, in situation where the CPU resources 
are limited. High accuracy of this approximation allows one to use it in a wide range of conditions, from Thomson to 
extreme Klein-Nishina limits. The approach adopted results also in some new analytic expressions representing known 
^_ ■ results in the Thomson limit. 
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1 Introduction 

rG • Since the first map of supernova remnant (SNR) RX J1713. 7-3946 in very high-energy gamma-rays ob- 

tained by H.E.S.S. (Aharonian et al. |2004|), there are few more SNRs with spatial distributions of VHE 
y-ray emission reported (see TeV Gamma-ray Source Catalog^ and reviews Rowell et al. (2005 1, Funk 
1200811 ). Surface brightness distribution of y-ray emission of astrophysical objects is an important possi- 
bility to test models of kinetics of astroparticles as well as dynamics of magnetic field and turbulence in 
astrophysical plasma. Thus, an actual and important task is the modeling of the respective surface brightness 
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distribution. 

Inverse Compton (IC) electron-photon interactions is one of the most important processes of gamma- 
ray production in SNRs. Application of the exact formalism of IC emission to the numerical calculation of 
SNR images requires quite powerfull CPU resources. In such a situation, an accurate approximation may 
Q\ . be of interest because it considerably reduces CPU time and also results in some new formulae in classical 

analysis (Sect. I2~2l and l2.31 l. 

[~^ . A common approach is to deal with IC emissivity for a given energy of the initial (monochromatic) 

photons (e.g. Jones [1968], Blumenthal & Gould [1970]) with assumption of a given shape of electron 
spectrum (power-law as a common choice). The resulting IC photon spectrum is then given by the integra- 
tion over the energy distribution of the field photons. Such an approach is essential for special cases of the 
field photon energy distributions. 

It is known however, that the black-body radiation field may be used for IC emission in many astrophys- 
ical objects. In particular, for SNRs under typical conditions, one may consider just black-body photons 
(either with single or with few different temperatures representing CMB/IR/optical radiation). In SNRs 
with no clear IR emission assosiated, the contribution from CMB photons dominates the role of infrared 
and optical photons (see discussion in Appendix in Lazendic et al. 1200411 ). The IR/optical photon fields 
may contribute typically 10%-15% of the ICfiux in such SNRs (Gaisser et al. Ifl998l . Baring et al. Ifl999ll ). 

In this note, we present an approximation for IC emissivity which may be applied to IC emission origi- 
nating from the black-body photon field with some temperature T. Since T is a parameter in our approach, 
the approximation may be used for calculation of IC radiation from different photon fields (CMB, IR, 
optical). The target radiation field around some SNRs (e.g. around Galatic center) may not be black-body 
and/or the contribution from IR/optical photons may dominate over CMB there (Porter et al. [2006], Hinton 
& Aharonian [ 2007 1). In cases when different components of the target radiation field may be approximated 
by a superposition of multiple Planck distributions with different T, our approximation may be used in a 
similar fashion. Namely, the overall IC emission will be the weigthed sum of single approximations, each 
with different value of the temperature. In cases when the initial radiation field may not be approximately 
described by a sum of black-body distributions, our approximation is not applicable. 

Another assumption in the present paper is the isotropy of the electron and photon fields. A thorough 
treatment of anisotropic IC scattering from cosmic-ray electrons is done by Moskalenko & Strong (|2000[). 

Our approximation is given in terms of an energy of incident electrons rather than in commonly used 
terms of the field photon energy. Such an approach opens the possibility for accurate modeling of IC emis- 
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sion of electrons with energy spectra being different from power-law. For example, if we are interested in 
electron energies around maximum possible values. It is known that contribution from electrons accelerated 
by the shock to E m - dx ~ 30 - 300 TeV is important in interpretation of the H.E.S.S. observations of shell 
SNRs. 

2 Approximation 

2.1 Overview of known formulae 

The spectral ditribution of the volume emissivity of (isotropic ally distributed) electrons due to IC process 
is (Jones IT968) , Blumenthal & Gould [19701 , Schlickeiser 1120021 ) 



P (E 7 ) = cE y I dyN(y) I den ph (£)o- KN (e 7 , e; y) 



(1) 



where y is the Lorenz factor of electron, N(y) is the spectral distribution of electrons, e and E y are the 
energies of photon before and after interaction respectively, n p h(e) is the (isotropic) initial photon energy 
distribution, 

o-kn (E y , e; y) = -^fi (q, rj) (2) 

is the angle-integrated IC cross-section, crj is the Thomson cross-section, 

G(q,rj) = 2q\nq + {\ + 2q)(\ -q) + 2r}q{\ - q), (3) 

E r r _ 4ey = eE y 

F(ym e c 2 — Ey) ' m e c 2 ' (m e c 2 ) 2 

The last term in G reflects the Klein-Nishina decline; recoil is embedded in the q parameter. The KN 

effect is more important than recoil for F > 1 (e.g. Baring et al. [1999]). Kinematic requirements result 

in (4y 2 y l < q < 1 (Blumenthal & Gould [1970|). Setting q to its minimum and maximum values limits 

energies of upscattered photons: 
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that simplifies to 
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in Thomson limit (T <c 1 ) and to 

E r , mi n = e (if T « 4y z ), E y , mSK = ym e c z ( 8 ) 

in extreme Klein-Nishina limit (F » 1). The condition q < 1 sets the minimum Lorentz factor 
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electron should have in order to scattter photon with energy e to energy E y . The function y m i a (E y ) may 
approximately be splitted in two parts 



E y /2 /(2e 1/2 ), for rj «C 1 
Ey/ (m e c 2 ) , for 77 » 1 



(10) 



switch from y^ m oc EJ to y m ; n oc E y is 



(The Klein-Nishina decline is negligible for 77 <K 1, Eq. [3]) The point where one could approximately 
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2.1.1 Method of approximation 

In some astrophysical environments, the initial photon energy field may well be represented by the isotropic 

black-body radiation 

1 e 2 

« P h(e) = 2 3 3 — — — - (12) 

n l Wc i exp(e/e c ) - 1 

with e c = kT. 

Let us re-write Eq. ([1]) in the form 
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dyN(y)p(y,E 7 ) (13) 



where the spectral distribution of IC radiation power of a "single" electron with Lorenz factor y is 

3cr T m 2 c 2 e c _ 2 2<? 4 e c _ 2 

p{y,E 7 )= y I(t] c ,t] q ) = -^7 iX»7c.J7o) (14) 



with the function I(rj c (E y ), rj (y, E y )) 

T(n„ n„\ = I ■ 

exp(j]/j] c )- 1 
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Let us introduce 

Gi(4)=2tfln2 + (l+24)(l-tf), (17) 

G2fo,J7o)=2J7o(l-«). (18) 

In the limit 77 — > 00, Gj, G2, G approach assymptotycally the values 

Gl,as = 1, G 2 , as = 277 , G as = 1 + 2rj . (19) 

The minimum value of 77, namely 77 m i n = j] , is given by the condition G(j] /t], rj ) - 0. Relation 77 m i„ = 770 
with definition of 77, Eq. (@J, yeild the formula for the minimum energy of the photon e m ; n which may be 
upscattered to the energy E 7 by the electron with Lorentz factor y: 

E y m e c 2 

(20) 
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In the limit 770 <k 77c that is equivalent to the Thomson limit T(e c ) <<c 1, the integral (fl~5T > may be found 
analytically 
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With decreasing of 77, G falls rather rapidly from G as to zero. Let G(t] /j], 77) be approximated by Heav- 
icide step function 

G(77 /77, 17) * (1 + 2t7 )//(t7 - 77J. (22) 



The integral I may then be approximately found as J(t7 c , ?7 ) ~ I~h(?7c, ?7o): 

exp (t7/t7 c ) - 1 



J H (77c, 170) = (i + 2770) f — rrV-r^ 

J exp(77/77 c )-l 

>7o k Z:5 J 

= (1 + 2t7 )t7 c — + Li 2 (exp(77 /77 c )) + 
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where the dilogarithm function 



1 
Li 2 (x) = - fj^-dt. (24) 
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Figure 1: The accuracy of the approximation (1321 as sum I] + Ii. Integrals i*i(l, x), J2O, x) (solid lines) 
are compared here with their approximations ( f30b and OTT ) (dotted lines). 



An accurate approximation of ( 1231 is 



2 / 9 

JhC^c^o) ~ -ri/c(l +2?7 )exp - — - 
6 \ 3r/ c 



(25) 



This expression restores J almost exactly for j] /t] c < 1CT 2 and r/ /?7 c > 1. However, In(r] c ,T] ) overesti- 
mates J(/7 C , 77 ) up to 2 times for 10~ 2 < rj /ri c < 1. This is because the deviation of G from Heavicide step 
function is important in this range. 

The original integral I has an important property. Namely, it may be splitted as I — I\ + 1 2 with 



I\(ji c ,t] )= — — — — flf/7 = 77c — — -dx, 

J exp(77/77 c )- 1 J expU)-l 
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(26) 



(27) 



= 77c?7o I rdx, 

where x - J] /j] c and G2* = G2/r] . As one can see, these integrals may be scaled: 
I\{arj c , arj ) = aJi(r/ c , ifo), Iiiarjc, ar] ) = a 2 I 2 (ri c , r/ ). 



(28) 



The possibility to scale these summands (integral J may not be universally scaled), is an important property 
which allows us to obtain an analytic approximation. It is also the reason of accuracy of approximation over 
the wide range of parameters. 

Really, using a - ?]~ l in the scale laws ( l28l i. one have 



-Ti07c,?7o) = rj c Ii(l,i] /i] c ), I 2 (r]c,ri ) = r} c I 2 (l,r} /r} c ). 



(29) 



7fc means that it is enough to check how accurate will be approximate expressions in approximation of 
just Ti(l, x) and l2(i,x), and we will know how accurate will be these approximations for any tj c and tj . 

Let us correct a bit exponential part in each summand of (f25l by introducing into exponent the second 
term of the form ci(?7 /^ c ) C2 where c\ and q are constant. These terms make the summands in ( EBT l to be 
more accurate in the representation of I\ and 1 \. It is important that the terms are also scaled in accordance 
to d28T i. With these corrections, we come to approximations 
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(30) 
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Figure 2: Accuracy of approximation d33l l. Integrals I ( fT5l l (solid lines), I\ and 1 2 d26l ). d27l > (dashed lines) 
and respective summands of approximation ( f33T > (dotted lines) versus 77 for a number of 77 c = 0.01, 1, 100 
(from below). 
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(31) 



The values of c\ and C2 are obtained by fitting the exact I\{\,x) and J2(l,x). They are compared with 
approximation on Fig. Q] which reveals good accuracy of obtained approximate formulae. 
The scaling property is a reason, why the sum of approximations < T30b and OTb 
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(32) 



is accurate to represent J in any regime, from Thomson to extreme Klein-Nishina. 

Eq. d32b is good for any rj c . We may suggest two a bit simple approximations for different ranges of tj c . 
If one are interested mostly in rj c < 100 (the case of IC emission of electrons accelerated by the forward 
shock in SNRs) then one can use an expression 



I(j)c, Vo) ~ -rrjc (1 + 3?7o) exp 
6 



2?7o 
3?7c 



4Wc 



1/2 



(33) 



Fig.|2]shows, that the second summand here (oc 3rj ) overestimates 1 2 in its power-law part. Nevertheless, 
this error is negligible for r] z < 100 (Fig. |2). If rj c > 10 is of interest, then one may neglect accuracy in 
exponential part of approximation of I\ and use the approximation (Fig. 01 



?(Jlc, Ho) ~ -rnc (1 + 2t7o) exp 
o 
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(34) 



Fig.|4]shows the accuracy of Eq.d33l in approximation of emission power ( fT4T > for electrons with ener- 
gies 10 12 -f 10 16 eV. 



2.2 'Delta-function' approximation 

Fig-0demonstrates that 'single' electron with Lorentz factor y - being scattered by all black-body photons 
- emits most of its IC radiation at photons with some characteristic energy E ym . 

One can introduce a 'delta-function approximation', in addition to the classical 'monochromatic ap- 
proximation' where electron is scattered by the photons with a fixed energy e (e.g. Schlickeiser [2002|). 
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Figure 3: Accuracy of approximation d34V Integrals I ( fT5l l (solid lines), I\ and It (1261 ). d27] > (dashed lines) 
and respective summands of approximation (f34b (dotted lines) versus 77 for a number of 77 c = 0.01, 1, 100 
(from below). 

Namely, one can assume that 'single' electron (scattered by all black-body photons) emits all of its IC power 
at photons with E ym : 

p(y, E y ) * p m (j)5(E y - E ym ) (35) 



where 



I' 



Pmiy) = I p{y,E y )dEy. 





(36) 



There are well known approximations for the total IC energy loss p m of electron in the Thomson (see ( |4TT ) 
below) and extreme Klein-Nishina limits (e.g. Sect. 4.2.3 in Schlickeiser [2002]). 
Our numerical calculations show (Fig.[4j that E ym may well be approximated by 



E 7 m(y) * E %imx (y, eb) 
where Ey m!LX {y, e) is given by ©. In the Thomson limit, this is 

Eymiy) ~ 4e c y 2 . 



(37) 



(38) 



Note, that in the classical 'monochromatic approximation', the average \E y ) - (4/3)e y 2 is used as an 
estimator for the energy of 'monochromatic' IC photons emitted by electron. 

2.3 Thomson limit 

The use of (f2Tb in ( TBI allows us to write down the expression for IC emissivity in the Thomson limit 
rjo <K 77 c . The spectral distribution of IC radiation power of a "single" electron with energy E = ym e c 2 is 



l F n _ °~T^ E y _ ne A el E y 
PT(T ' r) ~8/rVy 2 ~ 3hVE 2 ' 



tLy S ^y.lim 



(39) 



where E y ^ lm is a characteristic maximum energy defined below. This expression represents integration over 
all possible energies e of the seed black-body photons. 

The power pj(y, E y ) is increasing function of E y , while p(y, E y ) decreases rather rapidly after the 
maximum (Fig. |4j». Let's define an energy Zsy.iim by the condition 
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p T (y, Ey)dE y = p mT (y) 



(40) 
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Figure 4: The spectrum p{E y ) (TT4b calculated with integral I (solid lines) and with its approximation (l33l ) 
(dashed lines) for Lorentz factors y = 3 • (10 6 , 10 7 , 10 8 , 10 9 , 10 10 ) (from the left) that correspond to electron 
energies E — 1.5 • (10 12 , 10 13 , 10 14 , 10 15 , 10 16 ) eV. Crosses correspond to position of E JtmK (y, e c ). 



where 



Pmriy) = (4/3)ccr T wy 



(41) 



is the total energy loss of electron due to IC in the Thomson limit, u> — J en p \ l (e)de is the energy density of 
all initial photons. The definition j40l ) results in 
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= 2.65e c y 2 . 



(42) 



Note, that E y ^ m differs just a bit from E ym . dx (e c ) - 4e c y 2 . 



The volume emissivity of electrons distributed as 

N(y) = N y~ s for y min < y < y max , y min <s y tt 



(43) 



is 
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For £ y < 2? 7j *, the minimum Lorentz factor is y m ; n = £y /(2e* ), Eq. ( [Tol l, and we come to approximation 
of < n~3b in the Thomson limit 



P(E r ) = 



9.-2^. Jlls+D/2 
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(45) 



with the known slope W£y) « £ y . The value of e* may be fixed by comparison of ( 143b with e.g. 

expression (4.2.17) in Schlickeiser (|2002|). Namely, e* - J\{s)e c with 
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(46) 



Numerically, .71(1.8) = 0.665, 71(2) = 0.710, 71(2.2) = 0.755. 



3 Conclusions 

Numerical evaluation of the spatial distribution of the IC emission in SNRs requires essential computational 
resources because the IC volume emissivity (Q]) consists in two enclosed integrations (over initial photon 



and electron energies). We developed the approximation (132b of the spectral distribution of the IC emission 
power p{E y ) of electrons with Lorentz factor y which are interacting with the isotropical black-body photon 
field. Namely, Eq. d32l restores known results (Blumenthal & Gould [1970]) with high enough accuracy 
in any regime, from Thomson to extreme Klein-Nishina limits. It may be used for different astrophysical 

objects. For kTE y < 100(»i e c 2 J , i.e. in case of IC y-ray emission from electrons accelerated by the 
forward shock of SNRs, it is suitable to use a bit simple approximation ((33). In the Thomson limit, our 
approach results in a simple expression ( 1391 ). 

Our approximation may be used in situation when the initial radiation field may be approximated with 
the Planck function with some temperature T or when it may be represented by a superposition of the 
black-body distributions with different T. In addition, it assumes isotropy of electron and photon fields. 

The approximation is given in terms of an energy of incident electrons rather than in commonly used 
terms of the field photon energy (there are known approximation for the latter approach, see e.g. Schlick- 
eiser [2002]). Therefore, our approximation may be useful for analysis of the role of the electron spectrum 
with shapes different from power-law. 

The main idea behind our approach is the possibility to split initial integral into two parts which, contrary 
to the original integral, may be scaled. This scaling is the reason of accuracy of the approximation over the 
wide range of parameters, from Thomson to extreme Klein-Nishina regime. 

There is well known 'monochromatic approximation' where electron is scattered by the monochromatic 
photons with energy e (e.g. Schlickeiser [2002]). Fig.|4]shows that the spectral distribution of the radiation 
power of 'single' electron with Lorentz factor y scattered by photons distributed with the Planck function is 
peaked at some energy E ym (y). This allows us to introduce - similarly to the case of synchrotron emission 
- the 'delta-function' approximation for IC emission. In this approximation, Eq. d35l ). all radiated energy of 
electron is assumed to be at E ym (y). In the Thomson limit, E ym (y) » 4e c y 2 where e c = kT. In the classical 
'monochromatic approximation', the average \E y ) - (4-/3)e y 2 is used as an estimator for the energy of 
emitted IC photons. 

Our approach results in some new expressions which represent known results. Namely, Eq. ( |39l yeilds 
the spectral distribution of IC radiation power of a "single" electron and Eq. (l45l l represents the spectrum 
of IC emission from the power-law spectrum of electrons, in the Thomson limit. These expression account 
for integration over all possible energies of the seed black-body photons. They can be derived thanks to the 
possibility of analytical integration of J in the Thomson regime, Eq. d2"TV 
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